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SUB-FINSLER STRUCTURES FROM THE TIME-OPTIMAL CONTROL 
VIEWPOINT FOR SOME NILPOTENT DISTRIBUTIONS 


DAVIDE BARILARI, UGO BOSCAIN, ENRICO LE DONNE, AND MARIO SIGALOTTI 

Abstract. In this paper we study the sub-FinsIer geometry as a time-optimal control problem. In 
particular, we consider non-smooth and non-strictly convex sub-Finsler structures associated with 
the Heisenberg, Grushin, and Martinet distributions. Motivated by problems in geometric group 
theory, we characterize extremal curves, discuss their optimality, and calculate the metric spheres, 
proving their Euclidean rectiliability. 
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1. Introduction 


Sub-Finsler geometry is a natural generalization of Finsler geometry, sub-Riemannian geometry, 
and hence Riemannian geometry. In this paper, we introduce a very general notion of sub-Finsler 
structures: at each point of a manifold we consider a subspace of the tangent space endowed with 
a norm. Such a norm is not necessarily supposed to be strictly convex nor smooth, even away 
from the origin (the typical example is the norm). We will only assume that this norm changes 
smoothly with respect to the point of the manifold, in a suitable sense (see Section |1 . 1| ) . Particularly 
interesting examples are those norms that are “constant” with respect to the point, since these are 
the structures that appear in geometric group theory and the theory of isometrically homogeneous 
geodesic spaces. 

Lie groups equipped with sub-Finsler structures appear in geometric group theory as asymptotic 


cones of nilpotent finitely generated groups. Indeed, in [ Pan89 | Pansu established that, if we look 
at the Cayley graph of a finitely generated nilpotent group from afar, such a metric graph looks 
like a Lie group endowed with a certain left-invariant geodesic metric. Namely, the sequence of 
metric spaces {(Cay 5 (r), ^ps)}n€N of scaled down Cayley graphs of the nilpotent group P with 
generating set S and word metric ps converges in the Gromov-Hausdorff topology [pro99 1 towards 
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a Lie group that is stratified and nilpotent and is equipped with a certain explicit left-invariant sub- 
Finsler metric. We remark that such metrics come from structures that are never sub-Riemannian 
since the norms are characterized by convex hulls of finitely many points. 

Another setting where sub-Finsler structures appear is in the study of spaces that are isomet- 
rically homogeneous, i.e., metric spaces on which the group of isometries acts transitively. Us¬ 
ing the theory of locally compact groups and methods from Lipschitz analysis on metric spaces, 
[|Gle52| , BM46, MZ74, Ber89a, Ber89b], under the additional assumptions of being of finite dimen¬ 
sion, locally compact, and the distance being intrinsic, it has been proven that these spaces are 
sub-Finsler manifolds of the following type. Let G be a connected Lie group and H a compact 
subgroup. Let A C T{G/H) be a G-invariant bracket-generating subbundle of the tangent bundle 
of the manifold M := G/H. Consider a function F : A that is G-invariant and, for any p G M, 
F restricted to the vector space Ap is a norm, i.e., it is subadditive, absolutely homogeneous and 
vanishes only at 0. The sub-Finsler distance (also called Carnot-Caratheodory distance) associated 
with A and F is defined as 


d{p, q) = inf lengthy( 7 ), Vp, qeGjH, 


where the infimum is taken along all curves 7 tangent to A ioining p to o and for such curves 
length^( 7 ) ■= f F(i). 

Sub-Finsler structures also appear in different applications in control theory, as soon as one con¬ 
siders time-optimal driftless control problems where the controls enter linearly and satisfy polytopic 
constraints. As an example we can mention time-optimal control problems for three level quantum 


systems |BCC05 |. 

Our paper gives a contribution towards the understanding of the geometry of sub-Finsler spaces. 
Some natural problems are the regularity of spheres and of geodesics. For instance, it is an open 
question whether spheres, i.e., boundaries of metric balls, are rectifiable from a Euclidean view¬ 
point. Regarding geodesics, it is not even known if any pair of points can always be connected 
by a piecewise smooth length-minimizing curve. If this is the case, one would like to know if the 
number of such pieces is uniformly bounded. These are fundamental questions coming directly 
from the asymptotic study of nilpotent finitely generated groups. Indeed, there are conjectures 
about asymptotic expansions for the volume growth of balls of large radii that are related to the 
rectifiability of spheres and to the above-mentioned regularity of geodesics for the asymptotic cone, 
see pLD13| ]. 

The problem of finding length-minimizing curves in ‘constant-type’ sub-Finsler geometry can 
be locally reformulated as a minimum-time control problem for a system that is linear in the 
controls. Locally, on a manifold M, one considers k vector fields Xi,... defining the subspace 
of the tangent space and a symmetric convex body B C M^. The set B identihes the Finsler 
unit ball {uiXi + ... -|-| tt G B} in the distribution span{Ai,..., A^}. The problem of 
finding sub-Finsler shortest curves between two points pi,P 2 can be rewritten as the problem of 
finding an absolutely continuous curve 7 (-) : [0, T] —>■ M, together with a measurable function 
u{-) : [0,T] —> M*’, called control, that minimizes the time T > 0 and satisfies 


( 1 ) 


7(f) = tti(f)Ai(7(t)) H-h Ufc(f)Afc(7(t)), 

u{t) G B, 

7(0) =pi, 7(r)=p2- 


The purpose of this paper is to review some techniques of optimal control and to apply them to 
the study of some low-dimensional key examples: the Heisenberg, the Grushin, and the Martinet 
distributions endowed with norms whose balls are squares. 
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Each of these distributions is globally defined by the span of 2 vector fields Xi,X 2 : 


Heisenberg: 

M 

= 


II 

1 

X2 

— _i_ 

— '-'y ^ 2 ’ 

Grushin: 

M 

= m2. 


= dx, 

X2 

= Xdy, 

Martinet; 

M 

= M^, 


= dx + y^dz, 

X2 

= dy. 


The existence of time-minimizers for the problem (j^ in these three cases is a classical consequence 
of Filippov’s theorem. 

We consider the i°° and norm with respect to the above Xi,X 2 . More precisely, for every 
V G span{Xi(q'), X 2 (g)} C TqM, we consider the two quantities 


||u||oo := min{max(|r(;i|, |r(; 2 |) : v = wiXi{q) + W 2 X 2 {q)}, 

||u||i ;= min{|u;i| + |t(; 2 | : v = WiXi{q) + W2X2{q)}. 

Notice that for an arbitrary structure defined by vector fields Xi,X 2 the norm with respect to 
Xi,X 2 is the i°° norm with respect to + X 2 ) and \{Xi — X 2 ). In the case of Heisenberg 

group it is easy to see that the structure coincides with the structure in a new system of 
coordinates. For the cases of Grushin and Martinet we obtain, up to a multiplicative constant, the 
vector fields 

. , Grushin; M = M?, Xi = dx + xdy, X 2 = dx — xdy, 

Martinet: M = Xi = dx + dy + y'^dz, X 2 = dx — dy + y‘^dz- 

Regarding our first result, recall that a bang-bang trajectory is a finite concatenation of curves, 
called arcs, corresponding to a control that is constant with values in {(1,1), (1, —1), (—1,1), (—1, —1)}. 


Theorem 1. Consider the sub-Finsler structures for the Heisenberg, Grushin and Martinet dis¬ 
tribution defined by norm with respect to vector fields (|^) or (^). Then the length-minimizing 
trajectories are curves of two types: 

(i) one component of the control is constantly equal to I or —I, 

(ii) bang-bang trajectory. 

Moreover, the length-minimizing trajectories that are not of type (i) have at most 7 arcs. In addition, 
all curves of type (i) are length-minimizers. 


i) and type (ii) are not mutually exclusive. Moreover, it turns out that for 
(i) there exists a length-minimizing trajectory of type (ii) connecting the 


We remark that type 
every trajectory of type 

same two points. As a corollary, we deduce that any pair of points can be connected by an optimal 
bang-bang trajectory with at most 7 arcs. 

The proof is based on the classical Pontryagin Maximal Principle, for the description of extremal 
trajectories, i.e., trajectories that satisfy a first-order optimality condition. The bound on the 
number of arcs for optimal trajectories is obtained via second-order optimality conditions proposed 
by Agrachev and Gamkrelidze in [ AG9[)(| . 

In all the three cases we give a complete description of the sub-Finsler spheres. The case of the 
Heisenberg group was already studied in [ BLDlj ] with metric methods and it was proved that the 
sub-Finsler sphere is Euclidean rectifiable. Here we recover the shape of the sub-Finsler spheres 
and in addition we obtain the sub-Finsler front, i.e., the set of endpoints of extremals at a fixed 
time. 

We remark that apart from the Heisenberg distribution, which comes form a left-invariant struc¬ 
ture with respect to a group law, the other two examples that we study are not homogeneous struc¬ 
tures. Nonetheless, they come from projections of homogeneous structures on groups. Namely, the 
Grushin plane is a right/left quotient of the Heisenberg group equipped with a right/left-invariant 
strucure and the Martinet space is a right/left quotient the Engel group (which is the simplest 
stratified group of rank 2 and step 3) equipped with a right/left-invariant strucure. Gonsequently, 
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given a length-minimizer in Grushin (resp. Martinet) the curve in Heisenberg (resp. Engel) with the 
same control is length-minimizer as well. A feature of being quotients of nilpotent groups, with the 
induced projected structure, is that the distribution can be dehned by vector helds that generate 
a nilpotent Lie algebra. Nilpotency simplifies considerably the problem. For example, in our case 
bang-bang trajectories have piecewise polynomial coordinates. 

Thanks to the previous description of length-minimizing trajectories we are able to parametrize 
the spheres for the Grushin plane and the Martinet distribution, and we obtain as a consequence 
the following result. 

Theorem 2. Consider the sub-Finsler structures for the Heisenberg, Grushin and Martinet distri¬ 
bution defined by i°° and norm with respect to vector fields (|^ . Then the sub-Finsler spheres 
are Euclidean rectifiable, semianalityc and homeomorphic to Euclidean spheres. 


The semi-analyticity of spheres is interesting since it does not hold for the sub-Riemannian 
Martinet sphere, as proved in |ABCK97]. Indeed in that paper it is proved that the sub-Riemannian 
Martinet sphere is even not sub-analytic. 

We mention that there are a few other works that consider the view point of sub-Finsler geometry. 
A part from the previously mentioned ones, in the papers |CM06, CMW07| the authors study the 
sub-Finsler geometry, such as geodesics and rigid curves, in three-dimensional manifolds and in 
Engel-type manifolds. However, in those papers there is an assumption that is classical in Finsler 
geometry: the norm is assumed to be smooth outside the zero section and strongly convex. The 
present paper deals mainly with the case where these assumptions are not satished. Another notable 
paper is | GM13 l, in which the authors study the sub-Finsler geometry associated with the solutions 
of evolution equations given by hrst-order differential operators, providing one more setting where 
sub-Finlser geometry appears naturally. In a paper in preparation [ LDNG15 |, the authors study 
the Euclidean Lipschitz regularity of arbitrary left-invariant distances in a family of homogeneous 
groups, including the Heisenberg group. A signihcant remark is that the sphere with respect to the 
sub-Riemannian distance in the product of the Heisenberg group with the real line is a Lipschitz 
domain, while the sphere for the sub-Finsler structure on the same group admits a cusp, as 
observed in [ BLD13(] . 


1.1. Definitions. A function on M^, A: G N, is a norm if it is subadditive, absolutely homogeneous 
and vanishes only at 0. 

A sub-Einsler structure (trivialized and of constant-type norm) of rank at most A; on a smooth 
manifold M is a pair (/, || • ||) where || • || is a norm on and / 7>TMisa smooth 

morphism of bundles such that f{{p} x M^) C TpM, for all p € M. 

With every sub-Finsler structure we associate the distribution T) = f{E) and a norm on D 
defined by 

(4) ||n||sF = inf{||tc|| : f{p,w) = v}, for all n € Pp. 

A distribution V C TM is Lie bracket generating if (Lie(r(T>)))p = TpM, for all p G M. Here 
r(T>) is the collection of smooth vector fields tangent to V and, given a family H of vector fields, 
we denote by Lie (7^) and J-p the Lie algebra generated by F and the evaluation of the elements of 
F at a point p, respectively. 

Remark 3. The norm on M^, defined as above, in not necessarily Finsler in the classical sense, 
since it is not smooth away from the origin. As a consequence, even when F = TM the function 
V e->■ ||n||sF does not necessarily endow M with a Finsler structure in the classical sense. 

The notion of sub-Finsler structure introduced above, when some non-smoothness of the norm 
is allowed, can be seen as a particular case of the following more general class: A partially smooth 
sub-Finsler structure on M is a triple {E, || • ||e,/) where F is a vector bundle over M, || • Hf; is 
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a partially smooth Finsler structure on E (defined following Matveev and Troyanov [MT12]), and 
f : E ^ TM is a smooth morphism of bundles such that f{Ep) C TpM, for all p G M. The norm 
on the induced distribution can be defined in analogy with (®, replacing ||rt;|| by \\{p,w)\\e- 


Given a sub-Finsler structure with distribution V and norm || • || we say that an absolutely 
continuous curve 7 : [0, T] ^ M is horizontal if j{t) G and in this case its lenght is defined by 

^(7) = [ ll7(0llsFdb 
Jo 

We can then define the induced distance 


d{p,q) = inf{-^( 7 ) : 'y : [0,T] ^ M horizontal and 7 ( 0 ) =p,'y{T) = q}, 
which is well dehned and finite if the distribution E is Lie bracket generating. 


2. Sub-Finsler geodesics as minimizers of a time-optimal control problem 


Let M be a smooth manifold and (/, || • ||) a sub-Finsler structure on M. Notice that the bundle 
morphism / : M x —>■ TM determines k vector fields Xi,... ,Xk defined by Xi{p) = f{p,ei), 
where ei,... ,ek is an orthonormal basis for M.^. (Conversely, given any k vector fields Xi,..., Xk 
on M there exists a unique bundle morphism f : M x TM for which Xi{p) = f{p,ei).) 

The norm || - || identihes the set 

B := {w : ||tc|| < 1}, 

which is a closed, convex, centrally symmetric, and with the origin in its interior. (Conversely, any 
such a set is the closed unit ball of a norm on M^.) 

The problem of finding sub-Finsler geodesics, i.e., curves that minimize the length between 
two points p and g, can be reinterpreted as a time-optimal control problem, that is the problem 
of minimizing the time T > 0 for which there exist 7 : [0, T] —>■ M absolutely continuous and 
tt : [0, T] —>■ measurable such that 

r 7 (t) = ui{t)Xi{'y{t)) H-h Ukit)Xk{'y{t)), for almost every t G [0,T], 

(5) < u{t) G B, for almost every t G [0,T], 

[ 7 ( 0 ) =P, l{T) = q. 

Notice that the control function u might not be uniquely determined by the trajectory 7 , since 
the vector fields might not be linearly independent at every point. However, given the control u 
there exists a unique trajectory 7 satisfying j{t) = Yl^=i ^i( 0 ^j( 7 (^)) 7 ( 0 ) = p. 


2.1. Hamiltonian formalism and Pontryagin Maximum Principle. If a pair ( 7 ( ),n( )) is a 
time-minimizer for (^) , then it satisfies the first-order necessary conditions given by the Pontryagin 
Maximum Principle (PMP). Here we state a suitable version of the PMP for time-optimal control 
problem on a manifold M (see, for instance, |AS04, Corollary 12.12]). 

Define the Hamiltonian 


k 

( 6 ) n{X,p,u) := {\,f{p,u)) = ^Ui{X,Xi{p)), 

i=l 

for A G TpM, p G M, and u G For every u G M^, let u) be the vector field on T*M 

uniquely determined by the relation 

a{-,H{X,p,u)) = d(^x,p)'H{X,p,u), 

where a is the canonical symplectic form on T*M. 
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Define the maximized Hamiltonian 


( 7 ) 


if(A,p) := max{'H(A,p, tt) : u&B], 


Theorem 4 (PMP). Let be a time-minimizer for Problem (ID- Then there exist an 

absolutely continuous function A : [0, T] — >■ T*M and a constant Aq > 0 such that 

(i) X{t) € \ {0}, for every t e [0,r], 

(ii) the pair {X{t),'y{t)) satisfies the Hamiltonian equation 

(A(t), 7 (t)) = 'H{X{t),'y{t),u{t))), for almost every t G [0,T], 
which, in canonical coordinates, is 
dH dH. 

X{t) = - — {X{t),'y{t),u{t)), j{t) = —(A(t), 7 (t),u(t)), for almost every t G [0,r], 

(hi) H{X{t),^{t),u{t)) = H{X{t),^{t)) = Aq, for almost every t G [0,T]. 


If A('), 7 (-) satisfy for some u{-) and Aq the conditions (i), (ii), (hi) of Theorem^, we say that 
(A(-), 7 (')) is an extremal pair, that 7 (-) is an extremal trajectory, and that A(-) is an extremal lift 
of 7(-). 

For every vector field Y, if (A(-), 7 (')) is an extremal pair, then the function t {X{t),Y{'-f{t))) 
is absolutely continuous and its derivative satisfies 


( 8 ) 


L(Mt),y(y(t))} 


k 

{X{t),Y^ uj (t) [Xj , y] ( 7 (t))), 




for almost every t as it follows from the next classical computation. In canonical coordinates, 
thanks to point (ii) in Theorem ^ one has 


L(Mt),y(y(t)}} 


Y Clit)) 

( FlH \ ^ BY f)H 

-^(A(i), 7(i), u{t))j Y (7(t)) + X{tf—{j{t)) — {X{t),j{t),u{t)) 

F) f BY 

-Htf^il{t),u{t))Y (7(t)) + X{tf — {j{t))f{j{t),u{t)) 

{X{t),[f{;u{t)),Y]{^m- 


2.2. Second-order optimality conditions. Our aim is to recall necessary conditions for the 
optimality of an extremal trajectory whose corresponding control is piecewise constant. We refer 
to | AG90|] . (See also | AS03 , pig05(] .) 


Theorem 5. Let M be a smooth manifold and f : M x TM a sub-l°° structure on M. Let 

(j(-),u(-)) be an extremal pair for Problem (|5|) and let A(-) be an extremal lift of ^{■). Assume that 
A(-) is the unique extremal lift of 'y{-), up to multiplication by a positive scalar. Assume that there 
exist = To < Ti < T 2 <■■■< tk < tk+i = T and ,..., G such that u{-) is constantly 
equal to on (r^, Tj+i), for j = 0,..., K. 

Fix j = 1,..., K. For i = 0,..., K letYi = /(•, u*) G Vec(M) and define recursively the operators 
Pj = Pj-l = idvec(M); 

Pi = Pi-I O Vi G {j + 1,..., K}, 

Pi = Pi+i o e-(u+ 2 -u+i)ad(yi+i)^ Vi G {0,..., j - 2}. 

Define the vector fields 

Zi = Pi{Yi), ViG{0,...,i^}. 
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Let Q be the quadratic form 


(9) Q{a)= Y1 , 

0<i<l<K 

defined on the space 

{ K K 

a = (ao,.. .,aK) G = 0, ^aiZj( 7 (rj)) = 0 

i=0 i=0 

If Q is not negative semi-definite, i.e., if there exists a ^ W such that Q{a) > 0, then 7 (-) is not 
time-minimizing. 



3. Sub-£°° structures 

A choice of norm in that is of particular interest is the ^°°-norm, that is, 

||u^|| = |w^|oo := sup ItCjl- 

When the norm in the dehnition of sub-Finsler structure is the ^°°-norm we speak about sub-i°° 
structure. The corresponding time-optimal control problem (P) rewrites as 

( 7 (f) = ui{t)Xi{j{t)) H- \-Uk{t)Xk{j{t)), for almost every t € [0,r], 

(11) < |ui(f)| < 1, for every f G A:}, f G [0,T], 

[7(0)=p, 'y{T)=q. 

The maximized Hamiltonian (0) is 

(12) H{\,p) = |(A,Xi(p))| + ... + |(A, Afc(p))|. 

3.1. Switching functions, singular, abnormal, and regular arcs. With every extremal pair 
(A(-), 7 (')) and every j = 1,..., A: we associate the switching function 

t (^j(f) := {X{t),Xj{'y{t))). 

By formula (|^ we have that 

k 

(13) fij{t) = {\{t),Y, Uj{t)[Xi, Xj]{'y{t))) for almost every t. 

i=l 

The maximality condition (iii) of the BMP and ( |T^ imply that 

(14) |¥^i(A)|H -f l¥^fc(A)| = Ao, for all f 

and that, for all j = 1,..., A: and almost every t, 

(15) (pj{t) / 0 ^ Uj{t) = sign (pj{t). 

The restriction of an extremal pair (A('), 7 (-)) to some open nonempty interval I C [0, T] is called 

(i) an abnormal arc if ipj{t) = 0 on / for all j = 1,..., k; 

(ii) a If j-singular arc if (pj{t) = 0 on /; 

(iii) a regular arc if ipj{t) 0 for every f G I and for every j = 1,..., A:; 

(iv) a bang arc if the control n(-) associated with the trajectory is constant and takes values in 

{1,-iT. 
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Notice that a regular arc is a bang arc, but the converse is not true. Indeed, bang arcs can be 
singular (see Section . 

A bang-bang trajectory is a curve corresponding to a control that is piecewise constant with values 
{1,-1}^. In particular, a concatenation of regular arcs is a bang-bang trajectory, called regular 
bang-bang trajectory. When not specified otherwise arcs are assumed to be maximal, meaning that 
the restriction of the extremal pair to strictly larger open intervals is not an arc. 


Remark 6 . An arc is abnormal if and only if it is (/jj-singular for all j = 1,... ,k and if and only if 
Ao = 0. The latter equivalence follows from (14). In particular, if a trajectory contains an abnormal 
arc then the whole trajectory is an abnormal arc. 


4. Heisenberg group 


In this section we provide a description of the time-minimizing trajectories in the sub-£°° Heisen¬ 
berg group. These results have been previously obtained in | BLD13 | using methods of metric ge¬ 
ometry. The aim of this section is to illustrate how to exploit the geometric-control tools presented 
in the previous sections to recover such results. 

We consider the sub-.£°° structure on the Heisenberg group HI ~ determined by the vector 
fields 


(16) 


Xi = dx- 


X2 = dy + -d^ 


Let us introduce the vector field X^ = d^, which satisfies [Xi,X 2 ] = X^ and [Xi, X 3 ] = [X 2 , X 3 ] = 0. 
We use the notation from the previous section. Formula (|T^ gives immediately 

(17) (fl = -U2(P3, ^2 = UHP3, ip3 = 0, 


where </? 3 (t) := {X{t), X^l'ylt))). Notice that, since Xi,X 2 ,X 3 are linearly independent at every 
point, this is a reformulation in coordinates of the vertical part of the Hamiltonian system of the 
PMP. 

We characterize here below the abnormal, singular, and regular arcs for the associated time- 
optimal control problem. First, we show that there is no nontrivial abnormal trajectory. Next, 
we describe the structure of regular and singular arcs, showing that every nonconstant extremal 
trajectory is either a singular arc or a concatenation of regular arcs. Finally, we give a bound on 
the maximal number of regular arcs of a time-minimizer. 


4.1. Abnormal arcs. 


Lemma 7. The only abnormal arcs on HI are the constant curves. Consequently, no minimizer 
joining two distinct points is abnormal. 

Proof. From Remark we have = <p 2 it) = 0 for all t. By non-triviality of the covector A('), 
we deduce that </? 3 (t) 7 ^ 0 for every t. 

By the hrst two equations in (|T7|) , we get ui{t) = U 2 {t) = 0 for almost every t. Then the 
trajectory is constant and does not minimize the time. □ 

4.2. Singular arcs. 

Lemma 8. On HI the nonconstant trajectories that have singular arcs are exactly those for which 
there exists j G {1,2} such that uj is constantly equal to 1 or —1. All of them consist of a single 
singular arc and are time-minimizers. 

Proof. In what follows the roles of ui and U 2 are interchangeable. Consider an extremal trajectory 
that is not trivial and is <y 9 i-singular when restricted to an interval I, i.e., <^i = 0 on I. Because 
of Lemma 0, the trajectory does not have abnormal arcs, i.e., Aq 7 ^ 0. Hence, by (|TJ), ip 2 never 
vanishes on I. By (H), U 2 is constantly equal to 1 or —1 on I. From the first equation in ([T^) we 




have (/93 = 0 on /, and hence on the whole interval of definition of the trajectory. In particular, by 
dn we have that the whole trajectory is (/Ji-singular. 

Conversely, every trajectory corresponding to U 2 = ±1 constant and ui measurable with |ui| < 1 
has a (/?i-singular extremal lift with (p 2 = ^ and (pi = ips = 0 . 

Moreover, each of such curves 7 = (x, y, z) : [0, T] ^ H is time-minimizing since T = \y{0)—y(T) \ 
and \y\ < 1 for every trajectory of 7 = ^ 1 X 1 ( 7 ) -|- U 2 -^ 2 ( 7 )j with |tti|, |u 2 | < 1 - □ 

4.3. Regular arcs. 

Lemma 9. On H the trajectories that have a regular arc are regular bang-bang. Moreover, all arcs 
have the same length s except possibly the last and the first arc, whose lengths are less than or equal 
to s. At the junction between regular arcs the components ui and U 2 of the control switch sign 
alternately. 

Proof. Let I be an interval corresponding to a regular arc of the trajectory. Without loss of 
generality, ipi,ip 2 > 0 on I. Hence, by (|T^ we have u = (1,1) on I. Fix to G I- Two cases are 
possible: 

(a) (/? 3 (to) = 0: By ([ 1 ^ ) we have that ipi and ip 2 are constant along the entire trajectory, which 
is then a single regular arc. 

(b) y? 3 (to) 0: Denote by a the constant value of ip^. Using ( 0 ) we find 

ipi{t) = p)i{to) - a{t - to), P 2 {t) = q:) 2 {to) + a{t - to), \/t € I. 

Without loss of generality a > 0. Set ti = to + Pi{to)/cL- If the trajectory is defined up to time ti, 
then ifi and ip 2 are positive in the interval (to,ti). Also p>i{ti) = 0 . 

Since U 2 = 1 in a neighborhood of ti, we deduce that ipi is affine in a neighborhood of ti, with 
slope —a. Hence (pi < 0 < (^2 in a right-neighborhood of ti. Then ti is the starting time of another 
regular arc with control u = (— 1 , 1 ). 

Repeating this argument, backwards in time as well, we conclude that the extremal trajectory 
is the concatenation of regular arcs of length ip 2 {ti)/a = {ipi{to) + P 2 {tfi))/^ except possibly 
for the first and last arc, see Figure |l]. The switching occur alternately for ui and U 2 . □ 

The picture of the switching function, in the nontrivial case (b), is given in Figure |l|. 



Figure 1. The switching functions for the Heisenberg structure, when <^3 0. 


4.4. Bound on nnmber of optimal regular arcs. 

Proposition 10. A regular bang-bang trajectory with more than 5 arcs is not optimal. 
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Proof. Let us consider a trajectory with 6 bang arcs. By Lemma without loss of generality we 
can assume that the successive values of the control are 

(-1, -1), (-1,1), (1,1), (1, -1), (-1, -1), (-1,1). 

Denote the length of the internal bang arcs by s (recall that all arcs have the same, except possibly 
the first and last). 

We are going to apply Theorem ^ by taking j = 3. Let T3 be the third switching time. Since at 
T3 the function (p 2 switches sign, we have that ^ 2 (^ 3 ) = (A(r3),X2(7 (t 3))) = 0. 

Up to multiplication of A(-) by a positive scalar, we can normalize 993, which is constant, to —1. 
Hence, (fiirs) = s, which implies that A(-) is uniquely determined by the sequence of switching 
times. Set 

X+ = Xi+X2, X_=Xi-X 2 . 

Following the notations of Theorem we have 

Zo = ^ _ 2 sX 3 , 

Zi = = -X_ - 2 sX3, 

^2 = X+, 

^3 = AT-, 

Z 4 = = -X+ - 2sX3, 

^ gSad(X_)g^ad(-X+)('_^_^ ^ _ 2 gX 3 . 

A simple calculation shows that 

for = o-Qs = (Ti2 = CT23 = (T34 = (T45 = 2, 

<702 = <704 = <713 = <715 = <724 = <735 = 0, 

T'OS = 7'14 = <725 = “2. 

Decomposing the relation I]jLo“*^*(7(73)) = 0 on the basis {X+(7 (t 3)),X_(7 (t 3)),X3(7(r3))}, 
one gets 

— Q:0 + 0(2 — Q<4 = 0, —Qfi + 0:3 — (15 = 0, 2s( — Oo — Q<i — 0:4 — 0:5) = 0. 

Solving in ao,Q;i,a2) gives 

0:3 = —Q;2, Q ;4 = —OlQ + 0:2, Q<5 = —Qfi — 0 ( 2 - 

Notice that the relation a* = 0 is automatically satisfied. Then we can parameterize the 

space W appearing in the statement of Theorem I by ao,ai,a2, i.e., 

W = {(cio, 0 ( 1 , a 2 , — 0 ( 2 , —OLQ + 0:2 , —q:i — 0:2) | oxq, ai, 0:2 £ 

and write the quadratic form Q as 

Q[oiq, Qfi, 0(2) — 4q;oQ<i + 4q:oC12 — 4q:2. 

In particular, Q{1, 1,0) = 4 > 0, which implies that the trajectory is not optimal. □ 

4.5. Optimal trajectories and shape of the unit ball. Here we summarize the results obtained 
in the previous sections and we plot the unit ball in the Heisenberg group. 

Recall that once we characterize the controls ui{t) and U2{t) associated with an extremal trajec¬ 
tory, the trajectory itself can be recovered by solving the differential equation 

7(t) = 'Ui(t)Xi(7(t)) -b U 2 {t)X 2 {'y{t)). 
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(a) Singular arcs 



(b) Regular bang arcs 


Figure 2. Singular and regular arcs in the Heisenberg group 


By the coordinate expression 


(18) 


of the vector helds Xi, X 2 , this is equivalent to solving the system 

/ 

X = U\ 

^ y = U2 

= ^{U2X - Uiy) 


In particular, the trajectory is determined by its projection 7 onto the xy-plane, since the z coor¬ 
dinate of the trajectory can be found by integration. As it is well-known, it computes the signed 
area defined by the closed curve given by following 7 and then coming back to the origin along a 
line segment. 

As discussed in Lemma the singular trajectories correspond to the case when the control ni(t) 
is constantly equal to ±1 and U 2 {t) is free (or the symmetric situation). In Figure ^ we can see 
an example of such a curve when ui{t) = 1. Recall that these curves are optimal for all times and 
that, given one such trajectory, there exists a time-optimal bang-bang trajectory with at most 3 
bang arcs connecting the same endpoints. 

Regular bang-bang trajectories correspond to switching functions as in Figure |T|, where the 
controls switch sign alternately. These trajectories draw squares in the xy-plane as in Figure 

If such a trajectory has more than 5 bang arcs, then Proposition |l^ guarantees that the trajectory 
is not optimal. 

Notice that there exist time-minimizing curves of this kind with 5 regular bang arcs, as illustrated 
in Figure pb|. However, not all bang-bang trajectories with 5 bang arcs are time-minimizing. Indeed, 
if the underlying square is swept more than once, then the trajectory is no more a minimizer. 
Finally, let us also remark that for every minimizer with 5 regular bang arcs there exists a minimizer 
with 4 regular bang arcs joining the same endpoints (see again Figure [2bD. 

We stress that, by the classification of the previous sections, the only extremal trajectories 
connecting two distinct points on the same vertical line in the Heisenberg group are regular bang- 
bang. Once the shape of optimal trajectories is known, a picture of the Heisenberg sphere can be 


easily drawn. See Figure 3a and Figure 3b 


5. Grushin structures 

In this section we provide a description of the time-minimizing trajectories in two different sub-^°° 
structures in the Grushin plane. 

The classical sub-Riemannian structure on the Grushin plane is the metric structure on 
determined by the choice of the orthonormal vector fields 

X\ — X2 — xdy. 
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(a) Sphere (b) Front 


Figure 3. Sphere and front of the unit sphere in the Heisenberg group 


In other words the sub-Riemannian distance is characterized as follows 

r-T 

7 = uiXii'y) + U2{t)X2i'y), 7 ( 0 ) = po, 7 ( 1 ) =pi'> . 


d{pi,P 2 ) = hrf < / \Ju\ + u^dt 


KJo ' J 

The geodesic problem for this distance is equivalent to the time-optimal control problem defined 
by Xi, X 2 and u(t) G B, where B = {uf + <1} is the standard Euclidean ball. 


Due to the lack of symmetry of the Grushin structure, it is meaningful to consider two different 
sub-.£°° structures on the Grushin plane. 


5.1. The first structure. Consider the sub-^°° structures on determined by the vector fields 

( 19 ) Yi = dx + xdy, ¥2 = dx- xdy. 

Notice that Yi = Xi + X 2 and Y 2 = Xi — X 2 , so that we are considering the sub-£^ sub-Finsler 
structure associated with Xi, X2, up to a dilation factor. Similarly as in the Heisenberg group, let 
us introduce the vector field Y 3 = dy. 

The Lie algebra generated by Yi,Y2,Y3 actually satisfies the same commutator relations as in 
the Heisenberg group, namely 

(20) [Ti, Ta] = I3, [Hi, Ts] = [>"2, T3] = 0. 

The identities (pOf ) gives the same equations ( [T^ ) obtained in the Heisenberg case for the switching 
functions along an extremal trajectory 

( 21 ) ifi = -U 2 P 3 , ^2 = UiP3, P3 = 0 . 

In particular, (p 3 is constant. From the relation Y 2 — Yi = xYs we have the additional relation 

(22) = X(P 3 . 


In the case p 3 = 0, we get p 2 = Pi equals a constant, which is different from zero since the covector 
cannot be identically zero. Hence U\,U 2 are both 1 or both —1. In other words, every trajectory 
corresponding to = 0 is a horizontal line. Such curves are indeed time-minimizers. 

Let us then consider the case p 3 7 ^ 0. Under this assumption, if both pi and p 2 are vanishing then 
the trajectory is abnormal. In particular relation (22) implies that x{t) = 0 along the trajectory, 
hence we deduce that the trajectory is reduced to a point which is contained in the y-axis. 


Lemma 11. The only abnormal arcs on the sub-i^ structure on defined by the vector fields 
dD are the constant curves contained in the set {x = 0}. Consequently, no minimizer joining two 
distinct points is abnormal. 
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An analogous reasoning shows that there are no (/ 9 i-singular (resp. (/? 2 -singular) trajectories that 
are not abnormal. Indeed assume the trajectory is (/Ji-singular. Then (pi{t) = 0 for all t, that 
implies, by that U 2 = 0 for all t (recall that we are in the case ip^ 7 ^ 0). Thus ip 2 is necessarily 
identically zero and the trajectory is actually abnormal. The situation is analogue for (/ 92 -singular 
trajectories. 

Following the lines of Lemma one can then show that, if the trajectory is not abnormal, then 
it is bang-bang, and all internal arcs of a bang-bang trajectory have the same length s. 


Lemma 12. On the sub-£°° structure on defined by the veetor fields (IT) the trajectories that 
have a regular arc are regular bang-bang. Moreover, all arcs have the same length s except possibly 
the last and the first are, whose lengths are less than or equal to s. At the junetion between regular 
arcs the components ui and U 2 of the control switch sign alternately. 


5.1.1. Bound on number of optimal regular ares. Notice that if a bang-bang trajectory has an 
internal bang arc whose length is t, then ui and U 2 switch on the lines x = ±t/2 (see Figura Q). 

Regarding optimality, we claim that a bang-bang trajectory with 4 bang arcs is not optimal. 
Indeed, every extremal trajectory starting from y-axis is not optimal after it intersects again the 
vertical axis, as it follows by replacing the trajectory by its reflection along the y-axis. 


5.1.2. Optimal trajeetories and shape of the unit ball. The picture of the regular bang bang trajec¬ 
tories for this structure on the Grushin plane is given in Figure The corresponding picture of 
the unit ball is obtained in Figure 


y 


X 


Figure 4. 


Regular bang-bang trajectories for the Grushin structure (IT) 


5.2. The second structure. We consider now the sub-^°° structure on determined by the 
vector fields 

(23) Xi = dx, X 2 = xdy. 
and we introduce the vector field X^ = dy. 

The Lie algebra generated by Xi,X 2 ,X^ again satisfies the same commutator relations as in the 
Heisenberg group, namely 

[Xi, X2] = X3, [Xi, X3] = [X2, X3] = 0. 

Thus the identity (^ gives the analog equations ( [l7|) for the switching functions along an extremal 
trajectory 

(24) Pi = -U 2 P 3 , P2 = UiP3, ip3 = 0, 

In particular, p 3 is constant. From X 2 = XX 3 we have the additional relation 


P2 = XP3. 
13 











In the case 993 = 0, we get (/?2 = 0 and ipi equals a nonzero constant (otherwise the covector is 
identically zero). Reasoning as in Lemma we have immediately the following result 

Lemma 13. On the sub-i°° structure on defined by the vector fields ( |2^ the nonconstant 
trajectories that have singular arcs are exactly those for which ui is constantly equal to 1 or —1. 
All of them consist of a single singular arc and are time-minimizers. 

Let us then assume in what follows 0. 

Lemma 14. The only abnormal arcs on the sub-i°° structure on defined by the vector fields 
are the constant curves contained in the y-axis. Consequently, no minimizer joining two distinct 
points is abnormal. 

If 7 ^ 0 and the trajectory is not abnormal, then as in Lemma it is regular bang-bang, all 
arcs have the same length s except possibly the last and the first one, whose lengths are less than 
or equal to s. At the junction between bang arcs the components ui and U 2 of the control switch 
sign alternately. 

Moreover, on a regular bang-bang trajectory, U 2 switches on the line x = 0, since, if (^ 2 (i) = 0 
at a point t, then x{t)(p 3 = 0. Therefore if a bang-bang trajectory has an internal bang arc whose 
length is s, then ui switches on the lines x = ±s. Moreover, at ui-switching times the function ui 
goes from 1 to —1 if the switch occurs in the half-plane x > 0 while it goes from —1 to 1 in the 
half-plane x < 0 , since 

sign( 99 i) = -sign(u 2 </ 53 ) = -sign((^ 2 ¥? 3 ) = -sign(x(/?|) = -sign(x). 

5.2.1. Bound on number of optimal regular arcs. Regarding optimality, we prove the following 
lemma. 

Lemma 15. A regular bang-bang trajectory with more than 3 arcs is not optimal. If, moreover, 
the trajectory starts on the y-axis and it is optimal, then it has at most 2 arcs. 

Proof. First notice that, contrarily to what happens in the Heisenberg case, the role of the two 
vector fields Xi, X 2 is not symmetric. The replacement of {ui,U 2 ) by (—ui, —M 2 ) coupled with the 
reversion in the order of bangs, on the contrary, still yields a symmetric, equivalent, situation. This 
is a general fact, since it simply corresponds to reverse the parameterization of the curve. Looking 
at regular bang-bang trajectories (see Figure ^ one immediately recognizes that the proof of the 
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lemma can be given by looking at two types of bang-bang trajectories, whose successive values of 
the control are 

( 1 ,- 1 ), ( 1 , 1 ), (- 1 , 1 ), (- 1 ,- 1 ) and ( 1 , 1 ), (- 1 , 1 ), (- 1 ,- 1 ), ( 1 ,- 1 ), 

respectively. In the first case, one notices that reflecting the second and third bang arcs with respect 
to the y-axis yields another horizontal curve with the same length, which is not extremal. Hence 
the curve is not optimal. This argument also shows that regular bang-bang trajectories starting 
from the y-axis and with more than 2 bang arcs are not optimal. 

In the second case, let us apply Theorem ^ at the second switching time. One gets 

Zo = Xi + X2 + 2 sX3, Zi = -Xi+X2, Z2 = -Xi-X 2, Z3 = Xi - X2 + 2 sX3. 

Parameterizing the space W by the coordinates ao,ai we get that W = {(ao,ai, —cei, —oq) | 
ao,ai G M}. Normalizing 993 = 1 (uniqueness of the covector up to a positive factor is proved as in 
the case of the Heisenberg group), we write the quadratic form Q as 

Q(ao,ai) = 2q:o + 4aoai — 2ai. 

Since <5(1,0) is positive, the considered trajectory is not optimal. This concludes the proof of 
Lemma |T^. □ 

5.2.2. Optimal trajectories and shape of the unit ball. In this structure for the Grushin plane we 
have singular trajectories that are similar to the one obtained in the Heisenberg group, see Figure 
Let us stress that in this case the tangent vector of the curve is forced to be inside a cone 
whose width increases with the x coordinate. Regular bang bang trajectories from the origin are 

" y 

2 

.y = - 



Figure 6. 


Singular trajectories for the Grushin structure (^) 


illustrated in Figure |^. These trajectories lose optimality as soon as they reach the vertical axes. 
The picture of the unit ball in the Grushin plane with this structure is in Figure 

6. Martinet structures 


In this section we provide a description of the time-minimizing trajectories for two different 
sub-.£°° structures associated with the Martinet distribution. This is the easiest example where 
nontrivial abnormal minimizers appear. 

The classical sub-Riemannian structure on the Martinet space is the metric structure on 
determined by the choice of the orthonormal vector fields 

Xi =dx + X2 = dy. 


The sub-Riemannian distance is then 
d{pi,P 2 ) = inf < / \lu\ + uldt 


7 = uiXii'j) + U 2 {t)X 2 i'y), 7 ( 0 ) = po, 7 ( 1 ) = Pi 
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Figure 7. Bang-bang trajectories for the Grushin structure (|^) 




Figure 8. 


Sphere and front of the unit sphere for the Grushin structure (23) 


The geodesic problem for this distance is equivalent to the time-optimal control problem defined 
by Xi,X 2 and u(t) € B = {uf + u^< 1}. 

As in the case of the Grushin plane, due to the lack of symmetry, we are lead to consider two 
different sub-£°° structures. 


6.1. The first structure. Gonsider the sub-^°° structure on determined by the vector helds 

(25) Yi = dx + dy + y^dz, Y 2 = - dy + y^d^. 

Notice that Yj = Xi-\-X 2 and Y 2 = Xi — X 2 , so that we are considering the sub-.^^ Finsler structure 
defined by Xi and X 2 , up to a dilation factor. In analogy to the other cases, let us introduce the 
following vector fields defined by the commutators of the elements of the basis of the distribution 

(26) Ts := [>h, T 2 ] = Ayd ,, y4 := [Ti, [Ti, Ta]] = 49,, Ts := [^ 2 , [Ti, ^ 2 ]] = -49,. 

The switching functions associated with these vector fields and with an extremal pair (A(-), 7 (')) 
are 

i = 1,...,5. 

They satisfy the following system of differential equations 

(27) ifi = -'«2<y93, ip2 = UHP2., ^3 = + ^27^5, 

7^4 = 0 , (p5 = 0 . 


Remark 16. It follows from the bracket relations (^) that (/?4 and = — (/?4 are constants and 
we have = ycp^ = —yip^. In particular, if tp^ = 0, then (/J 3 is also constantly equal to zero, and 
</?i,(/j 2 are constant. 


Lemma 17. The nontrivial abnormal arcs on the sub-l°° structure on defined by the vector 
fields (p^ are the horizontal lines contained in the plane {y = 0}. 












Proof. Assume that the trajectory is not reduced to a point and it is abnormal on some interval I. 
In particular we have (pi{t) = = 0 for all t a I, while its control {ui{t),U 2 {t)) is not identically 

zero on I. From identities (p7|) one immediately gets that —U 2 {t)if^{t) = = 0. Hence, if 

we have that = y{t)ipi = 0 for every t (recall that 994 is constant), then y{t) = 0 for all t € /, 
otherwise (/94 = 0 and the covector is identically zero. In particular ui = ^2 on I and the trajectory 
is contained in a line {y = 0, z = zq}. □ 

Remark 18. Every costant trajectory (i.e., such that ui{t) = U 2 {t) = 0 on [0,T]) is also abnormal. 

The classification of extremal trajectories on Martinet is then reduced to regular and those that 
are singular with respect to exactly one control. 


6.1.1. Singular arcs. Let us now consider a singular arc. We show that in this case we can recover 
its (singular) control by differentiation of the adjoint equations. 

Indeed assume that the trajectory is (/ 9 i-singular, i.e., 991 = 0 on /, and we want to recover its 
associated control ui. Notice that \u 2 \ = 1 is constant and ipi = —U 2 ^z. By singularity assumptions 
= 0, that implies ipzit) = 0 for all t G /. We deduce that either (^94 = 0 or tti = U 2 on I. We 
have two possibilities:: 

(i) if <^4 = 0 then ui is free, 

(ii) if (^4 7 ^ 0 then ui = U 2 and the singular arc is also a bang arc, with no constraint on its 
length. Moreover y = 0 on such an arc. 

The situation with (/ 92 -singular arcs is perfectly symmetric. 


6.1.2. Regular arcs. Assume that both </9i(0), <^ 92 ( 0 ) 7^ 0 (the same holds for small times by conti¬ 
nuity). Because of Remark 16, we can assume that ^ 0 (otherwise the trajectory is made of a 
single bang arc). We want to show that 

(a) When (/9i(0)(/92(0) > 0 then the two switching functions are affine in a right-neighborhood 
of 0. 

(b) When </9i(0)<y92(0) < 0 the two switching functions are quadratic in a right-neighborhood of 

0 . 


On a bang arc the controls satisfy |ni| = |ti 2 | = 1 and thus we can differentiate the identity ( p7|) 
and get 


Cpi = -U2f>2. = -4u2(tti - ri2)(/94, ip2 = = 4tti(ui - tt2)(/94. 


In case (a) we have that tti = tt 2 = ±1, which implies (pi = f >2 = 0. In case (b) we have 
U 1 — U 2 = ±2 and consequently (pi and Lp 2 are constant and nonzero (recall that (/94 7 ^ 0 is constant). 
The equations for case (a) are 

pi{t) = (/9l(0) -h i</9l(0) = (/9l(0) - 'U2¥93(0)t, 

P2{t) = (/92(0) -h i</92(0) = (/92(0) 'Ul</93(0)t, 

= ¥^3(0). 

Notice that </ 93 ( 0 ) = ?/( 0)(/94 is zero if we start on the abnormal set. The equations for case (b) are 

Plit) = (/9l(0) - rt2(/93(0)t - U2{ui - U2)p>A — , 

ip 2 it) = (P2i0) + ril(/93(0)t Ul{ui - U2)(PA — , 

Psit) = ¥^ 3 ( 0 ) -¥ {ui - U2)(pAt. 

In particular, the constant (pA determines the convexity of the quadratic arc of the switching func¬ 
tions. 
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Lemma 19. A regular bang arc can enter in a singular arc only if the switching function is quadratic 
and has vanishing derivative at the switching point. 

Proof. Assume, for instance, that at some time to € I we have (pi{to) = 1 and ip 2 {to) = 0. Then 
the control ui{t) = sign(/?i(t) is constantly equal to 1 in a neighborhood Lto of and since (/?3 is 
continuous we deduce that ip 2 = uiPi, is also continuous in Utg. Since on the singular arc ip 2 = 0, 
we conclude. □ 

Next we discuss the possible behavior of the switching functions for regnlar arcs. Let us assume 
that (/Ji(0) > 0 and ^ 52 ( 0 ) < 0. In particular ipi and (/?2 are quadratic on a right-neighborhood of 0. 
We are reduced to three possible cases for the the switching function ipi\ 

- it never vanishes in the qnadratic part (we say that ipi is of type NI, for not intersecting), 

- it vanishes in the quadratic part and is tangent to the zero level [type T for tangent), 

- it vanishes in the qnadratic part and is transversal to the zero level {type I for intersecting). 
In Figure we picture the switching functions when cpi is of type NI, while Figures H and [Tl| 
correspond to type T and type I, respectively. 



Figure 9. 


Switching fnnctions for the Martinet structure (^) when ipi is of type NI. 



Figure 10. Switching functions for the Martinet structure (|^ ) when tpi is of type 
T. The relation between the length of the third and fourth bang arcs can be easily 
deduced from the expression of the switching functions. 

Assuming that there are only regular bang arcs along the trajectory (as it is always the case 
when ifi is of type NI or I) we have the following result. 

Proposition 20. The switching functions of a trajectory that has only regular bang arcs are peri¬ 
odic. 

The proof of Proposition is a simple consequence of the formulas of the switching fnnctions 
and Lemma |l^. When ipi is of type T, the only freedom is in the length of singular arcs. The 
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Figure 11. Switching functions for the Martinet structure (^) when ipi is of type I. 


order in which the switching occur is as in Figures in and 1^, up to the symmetry which sends 
(y?i, (/? 2 , ¥^ 3 , into (—(/?!, —^ 2 -, ¥^ 3 ) —^a) (which corresponds to a reflection y —>■ —y). 


Remark 21. It is easy to see from equations (^) and Remark that if (/?i is of type I then the 
y coordinate of the corresponding trajectory has constant sign. When restricting our attention to 
trajectories starting on the plane {y = 0}, we can then exclude that ipi is of type I. 


6.1.3. Bound on the number of regular arcs for optimal trajectories. The goal of this section is to 
prove the following result. 


Proposition 22. A bang-bang trajectory with at least one regular arc and with more than 7 arcs 
(either bang or singular) is not optimal. If, moreover, the trajectory starts on the plane {y = 0} 
and has more than 5 arcs, then it is not optimal. 

The proof works by applying several times Theorem |^. We should distinguish trajectories for 
which the switching functions are of one of the three types NI, T, and I. 

In order to reduce the number of cases to be studied, we use the fact that time-reversion and 
reflection y —>■ —y lead to trajectories with equivalent optimality properties. 

Switching functions of type NI. We start by considering tpi of the type NI, as in Figure |^. 

Lemma 23. A regular bang-bang trajectory of type NI with more than 5 arcs is not optimal. If, 
moreover, the trajectory starts on the plane {y = 0} and has more than 3 arcs, then it is not 
optimal. 


Proof. We prove the first part of the lemma by showing that concatenations of the type 
(28) (1,-1), (1,1), (-1,1), (1,1), (1,-1) 


are not optimal. All concatenations of 6 bang arcs, indeed, contain a concatenation of this type, 
up to symmetries (see Figure |l^. 

For concatenations of type (|2^), applying Theorem ^ at the second switching time, we get by 
computations as the one seen in the previous sections that the space W and the quadratic form Q 
in the statement of Theorem ^ can be written as 

W = {(aoj «!, 0, -ai, -ao) | ao,^ 1^}; Q{o(o,(^i) = 8{tial -\- t 2 aoai). 


Since Q is not negative semidehnite, the corresponding trajectory is not optimal. 

In order to conclude the proof of Lemma notice that, by metric considerations, if the trajec¬ 
tory starts from the plane {y = 0}, then it stops to be optimal at the middle of the third bang arc, 
see Figure 12. □ 















Switching functions of type T. We prove here the following result concerning trajectories corre¬ 
sponding to switching functions of the type T as in Figure 10. 


Lemma 24. A trajectory of type T with more than 7 arcs is not optimal. If, moreover, the trajectory 
starts on the plane {y = 0} and has more than 5 arcs, then it is not optimal. 

Proof. We hrst consider the situation where t^^^^ > 0 for every k. We notice that every concatena¬ 
tions of 8 arcs contains, up to symmetries, a concatenation of 6 arcs of the type 


(29) (1, -1), (1,1), (-1,1), (-1, -1), (-1,1), (1,1). 

(See Figure |T^. We are going to show that a concatenation as in (|^) is not optimal. 

For concatenations of type applying Theorem |5| at the third switching time (at which y = 0), 
we get that the space W and the quadratic form Q in the statement of Theorem ^ are written as 

W = {(ao, Oil, Oi2, —OiQ, ao — 02 , —ao ~ oii) \ ao, ai,a2 € M}, 

Q(ao,ai,a2) = 2(t2 — 2t %„^ a^ + 8t2aoai + 8tging®o«2 — ^t‘^mg'^2- 

Notice that Q is not negative semidefinite, since Q{e, 1/e, 0) = 2 e^{t 2 — 2tl^^^) + 8 t 2 > 0 for e small 
enough. Hence, the corresponding trajectory is not optimal. 

In the case where = 0, a concatenation as in (|2^) reduces to a concatenation of 4 bang arcs 

( 1 ,- 1 ), ( 1 , 1 ), (- 1 , 1 ), ( 1 , 1 ). 


Considering the following arc, we recover a concatenation as in (p 
tions as in the previous section show non-optimality. 


1), for which the same computa- 


The proof of Lemma 


starting from the plane {y = 0}, see Figure 1^ 


can be concluded as before by metric considerations for trajectories 

□ 


Switching functions of type I. We consider here trajectories corresponding to switching functions 
of type I as in Figure |ri|. Notice that such trajectories never cross the plane {y = 0}. 

We prove the following result. 

Lemma 25. A regular bang-bang trajectory of type I with more than 5 arcs is not optimal. 

Proof. By the same symmetry considerations as in the cases NI and T, we are left to prove that 
concatenations of the type 

(30) (1, -1), (-1, -1), (-1,1), (1,1), (1, -1), (-1, -1) 
and 

( 31 ) (- 1 , 1 ), ( 1 , 1 ), ( 1 , - 1 ), (- 1 , - 1 ), (- 1 , 1 ), ( 1 , 1 ) 

are not optimal (see Figure |l^. Notice than in both cases the trajectory is contained in {y < 0}. 

The application of Theorem ^ to the two cases is very similar leading (computing the quadratic 
form Q at the second switching time T 2 ) to the expressions 

Q{ao, 01 , 02 ) = -4(ti - tsjag - 4y(r2)aoai - 4(-2t2 + 2 t3 + y(r2))aoa2 - 4(2t2 - is - y{T2))al 

and 

Q{oq, 01 , 02 ) = -4(ti - tz)ol - 4y(r2)aoai - 4(-2ti -h 2t2 + y{T 2 ))ooa 2 - 4(ti - 2t2 - y{T 2 ))al 

respectively. In both cases, since y(r 2 ) < 0, one has that Q{e, 1/e, 0) = —4y(T2) -|- O(e^) is positive 
for e small enough. Theorem ^ then allows to conclude that the corresponding trajectories are not 
optimal. □ 
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6.1.4. Optimal trajectories and shape of the unit ball. Here we present the different pictures for the 
(x, y)-components of trajectories corresponding to switching functions of the form NI, T and I. The 
dashed lines correspond to the part of the trajectory which is no more optimal. 

In the case of trajectories of type NI we have the behavior in Figure [^. Trajectories of type T 


y 

- 2^ - 



X 


2ti 







Figure 12. Regular bang-bang trajectories of type NI for the Martinet structure (^). 

have singular arcs of arbitrary length (see two examples in Figure ^). Notice that the switching 
to singular always happens at points where y = 0, namely on the Martinet surface. The last case 



Figure 13. 


Regular/singular bang trajectories of type T for the Martinet structure (^). 


is given by trajectories of type I (see Figure ID- In this case the trajectory is contained in a strip 
Vo < y{t) < yi with eiter 0 < ?/o or yi > 0. In view of the optimality results one gets the following 


picture of the unit ball in the Martinet structure (25), see Figures 15 and 16 


6.2. The second structure. The second sub-Finsler Martinet structure on that we are going 
to consider is the sub-.£°° structure determined by the vector fields 

(32) X^ = d^ + y^d,, X2 = dy. 

We introduce the vector fields 


(33) X 3 := [Xi,X 2 ] =2yd,, X 4 := [Xi, [Xi, X 2 ]] = 0, X 5 := [X 2 , [Xi, X 2 ]] = 2d„ 
and the switching functions 

ipi{t) = (A(t),Xi(7(t))), 1 = 1,... ,5. 

The functions ipi satisfy the following system of differential equations 

(34) ipi = -U2<F3, 02 = UHP2,, 03 = UWi + 04 = 0, 05 = 0. 
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Figure 14. Regular bang-bang trajectories of type I for the Martinet structure (^). 



Figure 15. Unit sphere for the Martinet structure (^5|), view from the x-axis. 



Figure 16. Unit sphere for the 


Martinet structure (25), view from the y-axis. 


In this case the additional relations given by the bracket relations (^) are (p^ = yip^ and (^4 = 0. 
In particular system (^) reduces to 

(35) Pi = -U2Ps, P2 = UiP2., ^3 = mPb, ^5 = 0. 


Reasoning as in Lemma 17, we have the following characterization of abnormal arcs. 


Lemma 26. The nontrivial abnormal arcs on the sub-£°° structure on defined by the vector 
fields (^H) are the horizontal lines contained in the plane {y = 0}. 
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Indeed abnormal trajectories are described by the equations {y = 0,z = zq}. The fact that 
these trajectories are the same in the two Martinet structure under consideration reflects the fact 
that abnormal trajectories are independent of the choice of the frame (they depend only on the 
distribution). 

6.2.1. Singular arcs. The situation is not in this case symmetric with respect to ipi and (p 2 singular. 

Let us first consider a (/?i-singular arc, i.e., = 0 on I. Since U 2 = ±1 is constant and 0 = ifi = 

—U 2 if 3 , it follows that ¥? 3 (t) = 0 for all t € I. From (^3 = U 2 ^P 5 we deduce that = 0 and ui 
is arbitrary. In particular (p^ is identically equal to zero, which implies that the trajectory stays 
singular for all times. 

Let us then consider a (/ 92 -singular arc, i.e., (/92 = 0 on I. Since ui = ±1 is constant and 
0 = ip 2 = uiips it follows that (p 3 {t) = 0 for all t £ I. We deduce that 0 = (ps = U 2 V^ 5 - Hence, either 
(/95 = 0 with U 2 arbitrary, or M 2 = 0 on I with ip^ 7 ^ 0. In the first case the trajectory stays singular 
for all times, in the second case the singular arc is contained in the plane {y = 0 } and coincides 
with an abnormal trajectory on the interval I. 

6.2.2. Regular arcs. The analysis is similar to that of the previous section. We therefore omit the 
computations, which yield the following result. 

Proposition 27. A regular trajectory with more than 6 arcs (either bang or singular) is not optimal. 
The picture of the unit ball in the Martinet structure (|3^ is given in Figure 



Figure 17. Unit sphere for the Martinet structure (|3^ . 


7. Euclidean rectifiability and semi-analyticity of spheres 

By construction the i°° spheres for the Heisenberg, Grushin and Martinet structures studied 
above are homeomorphic to Euclidean spheres (5^ for Heisenberg and Martinet and for Grushin). 
Moreover these spheres are graphs of piecewise-polynomial functions. 

It follows that these spheres are Euclidean rectifiable and semi-analytic. Theorem |2| follows. 
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